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Abstract 

By using the convolution formula derived within the framework of relativistic quantum mechan- 
ics, we have examined the Fermi Motion effects on the ratios R p d/ PP = cr pd / (2a pp ) between the 
proton-proton (p-p) and proton-deuteron (p-d) Drell-Yan cross sections. We have found that in 
the small %i < 0.3 region, the Fermi Motion effect is less than 1% and our results for the ratios 
Rpd/pp = a pd /(2a pp ) agree well with the data at 800 GeV. In the large Bjorken X2 > about 0.4 
region, the p-d Drell-Yan cross sections can be influenced strongly by the Fermi motion effect. At 
120 GeV the predicted Fermi Motion effect can enhance the ratios Rpd/pp by about 20 % at X2 ~ 0.6 
and about a factor of 2.5 at X2 — > 1.0. Our results suggest that the Fermi Motion effect, along with 
other possible nuclear effects, must be included, in extracting the d/u ratios in the proton from 
the experiments on p-p and p-d Drell-Yan processes at large X2- 

PACS numbers: 



I. INTRODUCTION 



Since the asymmetry between the anti-up(u) quark and anti-down (d) quark distribu- 
tions in the proton was revealed by the New Muon Collaboration [l[ (NMC), a series of 
experiments jl|-[5[ on the di-muons (/i + /i~) production from the Drell-Yan[6| (DY) processes 
in p + p and p + d collisions had been performed at Fermi National Accelerator Labora- 
tory (Fermi Laboratory). The objective was to extract the d/u ratio of the parton distri- 
bution functions (PDF) in the proton. The information from these experiments and the 
measurements [l|, Q, 0] of deep inelastic scattering (DIS) of leptons from the nucleon have 
confirmed the NMC's finding, d/u > 1, in the region of low Bjorkin x < about 0.3. 

The ratio d/u > 1 signals the nonperturbative nature of the sea of the proton. Its 
dynamical origins have been investigated 0-23] rather extensively. Precise experimental 



determination of d/u for higher x > 0.3 is needed to distinguish more decisively these 
models and to develop a deeper understanding of the sea of the proton. This information 



will soon become available from a forthcoming experiment [241] at Fermi Laboratory. 

It is instructive to describe here how the p-d DY data were analyzed. In the analysis of 
the data of Ref . |E[ , the leading-order DY cross section of the p + N collision with N = p 
(proton), n (neutron) is written s 

da pN Atio 2 

= ^M/M(* 2 ) + flfrOflW] , (i) 

where the sum is over all quark flavors, a = 1/137, e q is the quark charge, }'n{x) is the 
parton distribution of parton q in hadron N, and M is the virtual photon or di-lepton mass. 
Here X\ and x 2 are the Bjorken-x of partons from the beam (p) and target (N), respectively. 
The DY cross section for p + d is simply taken as 

da pd da pp a pn 

+ TT^T- (2) 



dx\dx 2 dx\dx 2 dx\dx 2 

Clearly the effects due to the nucleon Fermi Motion in the deuteron are not included in 
Eq.([2]). The use of Eq.((2]) seems valid to a very large extent in the small x 2 region since the 
available data on nuclei indicate that DY cross section per nucleon is rather independent of 
the nuclear mass number A. However it is not clear whether it can be used to extract the 
d/u ratios in the proton from the upcoming experiments at large x 2 . 



In the DIS studies[30j, it is well recognized that the nuclear effects must be considered 
in extracting the parton distributions of the nucleon from the data. In particular it is 
mandatory to include nucleon Fermi motion effects by using various forms of convolution 
formula to express the DIS cross sections in terms of nucleon momentum distributions and 
parton distribution functions. In contrast, there exists very limited efforts to include the 



similar nuclear effects in analyzing the DY data on the deuteron [2|-|5| and nuclei [25H29 
As a step to improve the situation, it is necessary to investigate under what assumptions 
Eqs.([T])-(J2]) can be derived from a formulation within which the effects due to the internal 
motions of partons in the nucleon and nucleons in the deuteron can be defined rigorously 
and consistently. This is the main objective of this work. 

We will develop convolution formula for calculating the p-d DY cross sections within 
the relativistic quantum mechanics proposed by Dirac|31], as reviewed by Keister and 



Polyzou|32]. The same theoretical framework was taken in the studies of electron-deuteron 
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scattering 33, 34], electron- 3 if e scattering[35], and DIS on deuteron[36]. There exists con- 
volution formula for DY processes on nuclei , as given, for example, in Ref. 25] . However 



it is not clear how those formula can be related to the relativistic formulation considered 
in our approach. To explain clearly the content of our approach, we thus will give a rather 
elementary derivation of our formula with all approximations specified explicitly. We will 
apply our formula to analyze the available data at 800 GeV0] and make predictions for the 
forthcoming experiment 24J . 

In section II, we start with the general covariant form of the DY cross section and indicate 
the procedures needed to obtain the well known qq — > cross section a qq . The same 

procedures are used in section III to derive formula for calculating the p+N DY cross sections 
from a qq and the properly defined parton distribution functions of the nucleon. In section 
IV, we use the impulse approximation to derive the convolution formula for calculating p + d 
DY cross sections from a qq , f%, and the nucleon momentum distribution p p<1 {pn) of a fast 
moving deuteron with momentum p^. 

In section V, we describe the procedures for applying the developed formula to perform 
numerical calculations of p + p and p + d DY cross sections using the available parton 
distributions 37-41] and realistic deuteron wavefunctions j42|- IBI] . The results are presented 
in section VI. We will compare our results with the available data at 800 GeVj5[ and make 
predictions at 120 GeV for analyzing the forthcoming experiment [24] . A summary and 
discussions on necessary future improvements are given in section VII. 



II. COVARIANT FORMULA FOR DRELL-YAN CROSS SECTIONS 

The formula presented in this section are derived from using the Bjorken-Drell]46] con- 
ventions for the Dirac matrices and the field operators for Fermions and photons. To use 
the formula of Relativistic Quantum Mechanics given in Ref. [32], the plane- wave state \k > 
is normalized as < k\k >= 5(k — k) and the bound states |<& a > of composite particles, 
nucleons or nuclei, are normalized as < $ Q |$/3 >= 5 a ^. To simplify the presentation, spin 
indices are suppressed; i.e. \k a > represents \k a , X a > for a particle a with helicity A a . Thus 
the formula presented here are only for the spin averaged cross sections which are the focus 
of this paper. 

We consider the di-muons production from the DY processes of hadron (h) - hadron (T) 
collisions: 

h(p h ) + T{p T ) -)• fi + (k + ) + fT(kJ) + X h (p Xh ) + X T {p XT ) , (3) 

where Xh and Xt are the undetected fragments, and the four-momentum of each particle 
is given within the parenthesis. In terms of the partonic qq — > 7 — > [i + + //~ mechanism, 
illustrated in FigJTJ the covariant form of the di-muons production cross section can be 
written as 



(2tt) 4 1 dk+ dk- 1 

4[(p fc • Pt) 2 ~ mlml} 1 / 2 (2tt) 6 2E + 2E_ q~ 4 



do= AU _ _ r ^- 5 -^^ ^^±/^(fc f ,fc.)F^(p fc , Pr , g ) > (4) 



where and are the masses for h and T, respectively, E± = [k^+m 2 ^] 1 / 2 are the energies 
of muons y u ± , and q = k + + is the momentum of the virtual photon. The leptonic tensor 
is defined by 

r(k + ,kj = (2tt) 6 (2£ + )(2£_) < M_b'"(0)|0 >< Q\f(Q)\k + k_ > . (5) 
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Here the leptonic current is 

f(x) = e^(x)j^(x) , (6) 



where ipn(x) is the field operator for muon, and e = a/47to; with a = 1/137. By using the 
definitions Eqs.flH])-®, it is straightforward to get the following analytic form of the lepton 
tensor 

f^{k + ,k_) = -4e 2 [k1k u _ + k u + k» - ^(Jfe+.Jfe.+mJ)]. (7) 

Within the parton model, the hadronic tensor in Eq.(T4]) is determined by the current 
Jfi(x) carried by partons q or q 

F^(p h ,p T , q) = {2n) 6 {2E h )(2E T ) / dp Xh dp XT 5 A (p h + p T - p Xh - p XT - q) 

X h ,X T 

X < PhPT\Jfi(0)\p Xh dp XT >< Px T Px h \Ju(0)\p h p T > , (8) 

where 

M x ) = Z)[e g e]^ g (ar)y*V«(a?) • (9) 

Here ip q (x) is the field operator for a quark q with charge e q e; i.e e u — | and e-d = —\ for 
the up and down quarks, respectively. 

The above covariant expressions are convenient for deriving the formula which can express 
the hadron-hadron DY cross sections in terms of the elementary partonic qq — > cross 
sections. To get such formula, we first show how the elementary qq — > cross section 

can be derived from Eq.fllJ) with h — q and T — q. Explicitly, Eq.pj for the q(p g ) + q{pq) — > 
ji + {k + ) + n~{k-) process is 

- Wq . ff! miV „ i(k?^^r^ *->*gfa,». «)} • do) 
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The next step is replace the intermediate states \p Xh p XT > by the the vacuum state |0 > in 
evaluating the hadronic tensor Eq.([5]). We thus have 

*S(p*,Pff,ff) = (2nf(2E q )(2E q ) 

x < PgP g \M°)\° >< 0\J v (0)\p q ,Pq > 6*(p g +pg - q) . (11) 



Substituting parton current Eq.(J9J) into Eq. ffTTj) . the hadronic tensor Fjj£ then has a form 
which is the same as the leptonic tensor f^ u defined by Eqs.flSD-flSl) except that the momen- 
tum variables and the charges associated with the Fermion field operators are different. By 
appropriately changing the momentum variables in Eq.([7|), we obtain 

F$(p q ,p q , q) = -4[e,ef[^ + p£p£ - g^(p q ■ p- q + m 2 q )]5\p q + p q - q) . (12) 

By using Eqs.(j7J) and f[T2"j) . Eq. ffTUj) for the cross sections of q(p q ) + q{p q ) — >• /i + (k + ) +p,~(k_) 
can then be written as 

t no/ x (2tt) 4 1 dk+ dk- 1 _ 4 . N 

da = 4[(p q -p,y- m r* wr^E-™:7 5 {p « +p *- q) 

x8[k + ■ p q k„ ■ p q + k_ ■ Pq k + ■ Pq + m f + - 2 k ~ )2 + ml ^-P*) 2 ] . (13) 

It is convenient to express the q-q DY cross section in terms of the invariant function 
<? 2 = (Pq + Pq) 2 = (&+ + k ) 2 . After some derivations and accounting for the color degrees 
of freedom of quarks, we obtain 

da q KP q ,Pg) 47ra% 1 [q 2 - m^/4] 1 / 2 2 

ci^ = -7- e ^ [^-mg/4]V2 ^ - fc + ft) ) . (W) 

where N c is the number of colors. Taking N c = 3 and considering q 2 » m 2 and q 2 » m 2 , 
we then obtain the familiar form 

da qq ~{p q ,p q ) 4vra 2 „ 2 l 2 2 

^2 = "^rtg^ ~(P«+Pf) )• ( 15 ) 

The above expression is identical to the commonly used expression, as given, for example, 



in Ref. 41 



In the next two sections, we will derive formula expressing the p-N and p-d DY cross 
sections in terms of da qq (p q ,p q )/dq 2 given in Eq.(jl5p. To simplify the presentation, we 
only present formula for q in the projectile p and q in the target N or d. The term from 
interchanging q -h- q will be included only in the final expressions for calculations. 

III. p-N DY CROSS SECTIONS 

Eq.flU) for the p(p x ) + iV(p 2 ) -> /i + (*:+) + A*~(fc-) + A p (p Xp ) + Xat(p Xjv ) process is 
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where the hadronic tensor, defined by Eq.^SJ), is 

F^(P P ,PN,q) = (2n)%2E p )(2E N ){ £ J dp Xp dp x J 4 (p p + p N - p Xp - PXn - q) 

X p ,Xn 

x < PnP p \J»(0)\Px p Px n >< Px n Px p \Jv(0)\p p Pn >} ■ (17) 

Within the parton model, the DY cross sections are calculated from the matrix element 
< qq\Jfi(0)\0 >< 0| J u (0)\qq > which is part of the matrix element describing the annihilation 
of a q ( q) from the projectile p and a q (q) from the target N into a photon. To identify 
such matrix elements, we insert a complete set of qq states 



dPqdpq \PgPq >< p q p q \ 

between < p x p P \ and J M (0) (</„(0) and \ PpP N >) in Eq. (IT7j) and neglect any electromagnetic 
contribution from the undetected fragments X p and X?. We then have 



f J1v(PpiPNi q) = (2ir) 6 (2E p )(2E N ) J dp Xp dp XN 5 4 (p p + p N - p Xp - PXn - q) 

X p ,Xn 

X J d Pq dpq J dp'gdp'g < Pp \p qPXp >< P N \PqPX N >< PX p Pq\Pp >< PX N Pq\PN > 



as 



x < P q pM0)\0 >< 0M0)\KK > • 
By momentum conservation, the overlap functions in the above equation can be written 

< PqPx\p P > = < Px\bp q \ Pp > 

= <Pp P (Pq,Px)S(Pp - P q - Px) , (19) 

where b^ q is the annihilation operator of quark q. By using the definition Eq. ([T9l . Eq. (fI5j) 
can be written as 

F^(p P ,PN,q) = (2n) 6 (2E p )(2E N ) ]T J dp Xp dp x j\p p + p N - p Xp - PXn - q) 

Xp,Xw 

x [ J dp q \(j)p p {p q: p Xp )\ 2 8{p p -p q - p Xp )} 

X[J dpq\(j) pN {p q ,p XN )\ 2 5{p N - p q - Px N )} 

X < PqPq\M0)\0 >< 0\J„(0)\p q Pg > (20) 

The evaluation of the expression Eq. (l20l needs rather detailed information about the 
undetected fragments X p and X x because of the dependence of £ 4 (p p + p x —Px p ~Px N — q) 
on their energies p x and p Xn . To simplify the calculation, we follow the common practice 
to neglect the explicit dependence of the energies p Xp and p\ N of the undetected fragments. 
This amounts to neglecting the binding effects by setting p x ~ e\ and p\ N ~ e 2 to write 

^{P P + Pn - Px p - Px N ~ q) ~ S 3 (p p +p N - p Xp - p XN - q)8{p° p + p% - e x - e 2 - q°) , 

(21) 

where e\ and e 2 are constants. 
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We now define 

*ftfo) = E / dpxMpM^Px p )\ 2 5(p p - p g - p Xp ) (22) 
for the projectile p, and 

F p N (Pq) = E / d fe I ^ (Pff » P^iv ) I 2 $(Pn ~ Pg ~ Px N ) (23) 

x N 

for the target N. These two definitions and the approximation Eq. (l2~TT) allow us to cast 
Eq. f[2"Ul) into the following form 



F*?(Pp,PN,q) = (2n)%2E p )(2E N )Y, dp q dp q [Fl(p q )F« N (p,)} 

q J 

X[< PqPqM°)\ >< °\M°)\PqPq > 

xS 3 (p g + P g - q))6(p° p +p° N - ex - e 2 - q )} (24) 

We next make a reasonable approximation that the average energy e\ (e 2 ) in Eq. (|24p is 
the difference between the energy of the projectile (p) (target N) and the removed parton q 
( q); namely assuming 

5(p° p + p° N -e 1 -e 2 - q°) = 5((p° p - ei ) + (p° N - e 2 ) - q°) 

~ S(p° q + P°q ~ q°) ■ (25) 

Then Eq. (|24p can be written as 
F p ^(p P ,PN,q) = E/ dp q dpq[F« p (p q )Fl(p q )]^ 



x{(2n) (i (2E q )(2E q ) < p q p q \J p (0)\Q >< 0| U0)\p q Pq > 5 4 (p q + p q - q)} . 

(26) 

The quantity within the bracket {...} in the above equation is just the hadronic tensor 
Fjj£(p q ,p q ), defined in Eq.( !TT|) . for the qq system. We thus have 



Fff(Pp,PN,Q) = E/ dPqdMFl(p g )F^ N (pq)]^F^ Pq ,p q ) . (27) 



Note that the above simple expression is due to the use of the approximations Eqs. (12 ip and 
( 12 5p . The binding effects on the partons in the nucleon and the undetected fragments X p 
and Xn are not treated rigorously. If we depart from these two simplifications, we then need 
the spectral function of the nucleon in terms of parton degrees of freedom to calculate DY 
cross sections. Such information is not available at the present time. 
Substitute Eq. (l27]) into Eq. ffTB"]) . we then have 

1 rlk rlk 1 

X {(2^^ r " ( * + '*-« (P " P5)} < 28) 



The quantity in the bracket {..} of the above equation is precisely what is in the bracket 
{..} of Eq. tflOl) for the qq — > process. Accounting for the difference in flux factors and 

extending Eq. (f28|) to include the q -H- q interchange term, the full expression of the p-N DY 
process is 



dq 2 



€gg -Pq) 2 ~ m*} 1/2 E P E N d^(p q ,p q ) . 
4[(p P ■ Pn) 2 — mhn 2 ^)] 1 / 2 E q E q dq 2 



(29) 



where — is the q-q DY cross section, as defined by Eq.flTSj). 

We now examine the physical meaning of the functions Fjj (p g ) and Fj> N (pq), defined in 
Eqs. fl22|) - f[23|) . The probability of finding a quark q with momentum p q in a nucleon state 
\p p > is defined by 

Inserting a complete set of states 1 = / dp Xp \px p >< Px p \ into the above equation and using 
the definition Eq. (fT9~l) . we then have 

,^ J2x p fdp Xp < P P \bl q \px v >< Px p \b? g \p P \ > 

<m> 

E Xp I dp q | 0p- p {Pq, PX p )\ 2 S (p p -Pq-pX p )S {p P ~ Pq ~ PX p ) 

With our normalization < p p \p p >= 5(p p — p p ), the above equation becomes 

= 1 igFiu 

= 53 / d Pq\^Pp(.P^PX P )\ 2 KPp-Pq-PX p ) (31) 

Clearly Fi p (p q ) of Eq. (f22l) is identical to Pf p (p q ) given above and is the probability of finding 
a quark with momentum p g in a nucleon moving with a momentum p p . Thus Fi p (p q ) is the 
parton distribution function in the nucleon. 

IV. p-d DY CROSS SECTION 

For the p(p p ) + d(p<j) — >■ + DY process, Eq.fllJ) gives 

dffPJ = w^i^^ l A^^ r{k+ - k - )F » i[p " pd < q)] • (32) 

We assume that the p-d DY process occurs only on each of the nucleons in the deuteron, as 
illustrated in FigfJJ In this impulse approximation, the hadronic tensor for a deuteron target 
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FIG. 2: Impulse approximation of p + d Drell-Yan process. 



can be obtained by simply extending Eq. (1171) for the p+N to include a spectator nucleon 
state \p s > in the sum over the final hadronic states. We thus have 

Fj£(Pp,Pd,q) = {2nf{2E p )(2E d ) J2 J d Psdpx p dpxJ A (p P + Pd-px p -px N -Ps- q) 

X p ,Xn 

x[< ® p J>p\Jp(G)\px,PX N Pa >< PsPX N Px p \Ju{0)\Pp®p d >} , (33) 

where |$ Pd > is the state of a deuteron moving with a momentum p d . Within the framework 
of relativistic quantum mechanics 32], we can expand the deuteron state in terms of two- 
nucleon plane- wave state \pxP2 > 

\® Pd >= J dp N dp 2 ® Pd (p N )5(pd ~Pn - P2)\pnP2 > • (34) 

Keeping only the contributions which are due to a parton in \p^ > of the above expansion 
and a parton from projectile state \p p >, we have 

< % d P P M°)\Px P Px N Ps > 

= J dp 2 J dp N $*{p N )5{p d -p N -p 2 ) < PPn\Jh{0)\px p Px n >< P2\p s > 

dp N §*(p N )5(p d -p N - Ps) < PPn\Jp(0)\px p Px n > • (35) 

By using Eq.(l35l). Eq.( l33|) can then be written as 

(2E p )(2E d ) 



F^(p P ,Pd,q) = J dp N \% d (p N )\ 



(2E P )(2E N ) 

x{(2it) 6 (2E p )(2E n ) x ^ | dp Xp dpxJ 4 (p P + Pn - Px p - Px N - q) 

X p ,Xn 

x[< PnP p \J^)\Px p Px n >< Px n Px p \Ju(0)\p p Pn >}} • (36) 

We see that the quantity within the bracket {..} in the above equation is identical to 
FP^(p p ,p N , q) of Eq.flPTj). We then have 



F^(p P ,Pd,q) = J dp N p Pd (p N ) 



E„E, 



P d F^(p p ,p N ,q) 



E„E 



N 



(37) 
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where 



p Pd (p N ) = \% d (p N )\ 2 . (38) 

By using Eq. (l3"I|) . one can show that p Pd ij>N) is the nucleon momentum distribution in a 
deuteron moving with a momentum p d : 

Pva (PN) = <$ ft • (39) 

We will present formula for calculating p Pd (j>N) in section V. 
By using Eq. (l3"7|) . Eq. (!32|) becomes 

( 27r ) 4 V- / ^ „ ^ N ^ 



daPd = Mv-vrf-rnhnM* JdPNP Pd (PN) 



4[(p ■ P„) 2 - m> 2 // 2 ™ — ^ 

X ^fefe^^)^^^>- (40) 

The quantity within the bracket {..} of the above equation is exactly what is in the bracket 
{..} of Eq.( jT6l) . Accounting for the difference in flux factor, we obviously can write 

d pd (p, Pd ) ^ f„ ^ , j(p ■ Pn) 2 ~ m 2 p m 2 N } 1/2 E p E d d? N {p,p N ) 

~ tkj PNPPi[PN \{p-p d ?-™> 2 d )Y /2 E p E N ' (41) 

where dP j^f^ is given in Eq.f l29|) . Substituting Eq.f l29|) into Eq.t HI]) . we have 

d pd (p, Pd ) ^ u - f [(p, ■ Pq? ~ ™t} 1/2 E p E d 



J d PNP Pd {PN){Yl J dPqdPq 



dq 2 ~ N ±;,J ™W'^J Fq [{p.p d f-mlm 2 d )]^E q E, 

dq 



* K(p q )F« N (p,) + F* N {p q )FHp q )][ rri } • (42) 



Clearly the input to our calculations of the p-d DY cross sections are the nucleon momentum 
distribution p Pd (pN), the parton distribution functions Fi{p q ) and F§(p q ), and the elemen- 
tary partonic qq — > p + p~ cross section — ^pr^ defined by Eq. f[T5"j) . 

Here we note that within the impulse approximation, as illustrated in FigfSJ the binding 
on the nucleon in the deuteron is treated accurately if p Pd (pW) is calculated from realistic 
nucleon-nucleon potentials which give high precision fits to the NN scattering data and the 
deuteron properties. In the considered framework of Relativistic Quantum Mechanics, the 
nucleon emitting a parton is on its mass-shell, while the interacting two-nucleons system 



is off the energy-shell of the deuteron: \Jp± + m 2 N + yp 2 2 + m ?v \JPd + m l- Precisely 
because of this choice of formulation, the parton distribution in Eq.( l42l) is the same as that 
in the free nucleon. If one chooses a theoretical framework within which the nucleon in the 
deuteron can be off its mass-shell p 2 ^ ttl 2 n , then this simple relation is lost and one needs 
addtional theoretical assumptions to extract the d/u ratios in the proton from p-p and p-d 
DY data. 
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V. NUMERICAL PROCEDURES 

In this section, we develop numerical procedures for applying the formula presented in 
previous sections to investigate the Fermi motion effect on the ratio R p d/ PP = cr pd / (2a pp ) 
between the p-d and p-p DY cross sections. Our first task is to relate our momentum 
variables p p , px and q to the variables used in the analysis|5j of the available data. This will 
be given in the first subsection V.A. The procedures for calculating DY cross sections are 
given for p-p in subsection V.B and p-d in subsection V.C. The calculations of momentum 
distributions P Pc1 {j>n) of a moving deuteron are described in subsection V.D. 



A. Kinematic variables for DY cross sections 

It is common jEf to use the co-linear approximation to define the parton momentum: 

Pq P = xiPp , (43) 
P qT = x 2 p T , (44) 

where p qp (p qT ) is the momentum of a parton in the projectile (target), and X\ and x 2 are 
scalar numbers. The momentum q of the virtual photon in the qq — > 7 — > / u + /i~ can then 
be written as 

q = Pq+ Pq 

= x x p v + x 2 p T . (45) 

In the considered very high energy region E p > 100 GeV, the masses of projectile (pp 1 = m£) 
and target (p? = m?) can be neglected and we have 

s = {p p + Pt) 2 ~2p p -p T , 
p T -q~ XiPp ■ p T , 
P v - x 2 p p ■ Pt ■ 



The above relations give 



xi , (46) 



s 



x 2 ~ • (47) 

s 

It is most convenient to perform calculations in terms of x\ and x 2 in the center of mass 
(cm.) system in which the projectile is in the z direction and the target in -z direction: 



Pp = UP 2 + m 2 p , 0, 0,p) ~ (p, 0, 0,p) , (48) 



Pt = (yp 2 + m 2 T , 0, 0, -p) ~ (p, 0, 0, -p) . (49) 

With the choices Eqs.pSJ-PSJ, we 

s = (p p +Pt) 2 ~ 4p 2 , 
M 2 = q 2 = (xiPp + x 2 Pt) 2 ~ ^xix 2 p 2 . 
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The above two equations leads the simple relation 



M 2 

%ix 2 ~ • (50) 

s 

By using Eqs.(l46l)-(l471) and Eqs.( HH|) -(H9l). we can define a useful variable xf 

2(Pt-Pp) ■ q 



Xf — Xl — X2 



s 

■ q 2y/sp p ■ q 



s 

In the notation of Ref. jH[, we write 



(51) 



7 

XF ~ ' (52) 

where pj = p p ■ q is clearly the longitudinal momentum of the intermediate photon with 
respect to the projectile in the cm. frame. Experimentally, s, M, xp, and da j [dMdxp) are 
measured |47|]. With the relation Eq. (1501) . we certainly can determine the corresponding x±, 
X2 and da/(dxidx 2 ). We thus will only give the expression of da / (dxidx 2 ) in the following 
subsections. 

Before we proceed further, it is necessary to discuss how our formula can be expressed 
in terms of X\ and x^- As given in Eqs. fl43l) - fj4ll) . these variables define the fraction of 
the hadron momentum carried by the emitted parton. For the projectile proton, this is 
unambiguous X\ = p q /p p = p\jp z v since we choose p in z-direction. For the deuteron target, 
X2 could be identified with the deuteron or the nucleons inside the deuteron since in the 
above derivations hadron masses are all neglected. Since the impulse approximation is used 
in our derivations, we identify x 2 with the mean nucleon momentum < Pn > in the deuteron. 
As will be seen from the realistic nucleon momentum distribution to be presented in Figj3j 
< Pn >~ |Pd|/2 as expected naively. 



B. Calculation of p-p DY cross sections da pp / dx\dx 2 

We now note that with the variables X\ and x 2 defined above the flux factor associated 
with Eq. fl29l become 1. Substituting Eq. ffl5|) into Eq. ff29l . the DY cross section for p{p p ) + 
p(Pn) — > is then calculated from 

x ^ e2 Ai 2 -(Pi + p^))- ( 53 ) 

In the chosen center of mass frame, defined by Eqs. (|48l) - (l49l) . let us consider q in the target 
proton moving with a momentum p^ = (p z N , 0, 0, p z N ). In the precise co-linear approximation, 
only the z-component of the q momentum is defined by p z N . We thus write Pg = (pq,T,pf) 
where 

P\ = %2P Z N , (54) 
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and the transverse component pq t T can be arbitrary. The integration over the q momentum 
distribution in the target iV can then be written as 

/ dp q F} N {p q ) = J dpi J dp q , T F«{pl,p qtT ) 

= PN J dx 2 J dpg, T F^(x 2 p Z N ,Pq,T) 

= J dx 2 f* N (x 2 ) , (55) 

with 

fL fa) =P Z N J dp q , T F* N (x 2 p z N , Pg iT ) . (56) 

We identify fp N (x 2 ) with the parton distribution functions (PDF) determined by several 
groups [37H41|. To compare with the results of Ref. 0], we use the PDF of CETEQ5mj4l|. 
Similarly, we can define for the projectile proton 

jdp q Fl{p q ) = Jdx 1 f^x 1 ), (57) 

where (xx) is defined by the same Eq. (l56l) with q replaced by q and x 2 by x±. By using 
Eqs.( l56|) and f l57|) . Eq.(l53l) can be written as 

dPN( ^ PN) = E/ dx^lf^f^ + f^f^fa)] 

x^fe 2 q S(q 2 -(p q + pq) 2 ). (58) 

By integrating over dq 2 , the above equation leads to 



which is the same as Eq.flT]) used in the analysis of Ref.pj since (p q + Pq) 2 = q 2 = M 2 for 
the partonic process qq —> 7. 

We only consider up (u) and down (d) quarks in the proton. Eq. (1591) for the p-p DY cross 
section then obviously takes the following form 

d^(p, Pp ) _ W r 4 (/ « (a;i)/ « (a;2) + / « (xi)/ « (a;2)) 



d Xl dx 2 9M 2L 9 WpV JJpp 
1 



+lif d M)ftM) + f d P faKfa))\ • ( 6 °) 



C. Calculation of pd DY cross sections of da pd / dx\dx 2 



From Eq. (|42p . we see that the calculations of p-d DY cross sections involve the convolution 
of the parton distribution associated with a nucleon with momentum p^ over the nucleon 
momentum distribution p Pc1 {pn) in a moving deuteron. To see this more clearly, we consider 
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the contribution from an anti-quark in a nucleon N of the deuteron target and a quark in 
the projectile proton. By using Eq. (fT5|) . this particular contribution to Eq. (|42|) is 

d pd (p, Pd ) f f f l(p q -PtY-mt] 1/2 E N (p)E d (p d ) ^ 
dq 2 J Pq J Pqi [(p-p d ) 2 -m 2 p m 2 )Y/ 2 E q (p q )E q (p q )^ 

f - Airct 2 

x J dp N p Pd (p N ) [f*@ q )f* N (Pq] -^- e l 6 (V 2 ~ (Pa + Pa)) • ( 61 ) 

Bu using the variables X\ and x<i defined in subsection V.A, the ratio between flux factors 
in the {..} bracket of the above equation is close to 1. We thus only need to consider 



d (p,Pd) 
dq 2 



q j u r ,, i dp N p Pd (p N )[f*(p q )f« N (p. 

x^e^(g 2 -(p 9 + p 5 ) 2 ) (62) 

We now define the momentum fraction rc^ of a parton associated with a nucleon with 
momentum Pat in the deuteron 

Pq = %2P Z N- ( 63 ) 

By using the definitions Eq.(l63l) for and Eqs.(l56 l) -( !57l) for parton distributions, Eq.(!62l 
can be written as 



d pd (p, Pd ) 
dq 2 



dxx / dx% / dp N p Pd (p N )[f«(xi)f« N (x%)) 



x ^ e ?^ 2 -(^ + ^-) 2 ) ( 64 ) 

Similar to the pp case, p d = {p%Pd,T = 0) is chosen to be on the z-direction of the 
cm. frame. As discussed above, we need to find a way to relate to 22 which is determined 
by experimental variables M, s and Xf through the relations: X1X2 = M 2 / s and Xf = 
X1—X2 = p7J(-\/s/2). In the impulse approximation we have used to derive the p-d DY cross 
section, it is assumed that the parton is emitted from one of the nucleons in the deuteron, 
as illustrated in Figj2j Thus it is reasonable to choose p q = X2p z ave where p ave is the averaged 
nucleon momentum in deuteron. Combining this assumption and Eq.( |63l) . we then have the 
following relation 

x n = X2 P^ . (65) 
Pn 

In section VI, we will determine p ave from investigating p Pc1 (pn) generated from realistic 
deuteron wavefunctions. 

By using Eq. flo^]) to change the integration variable, we can write Eq.flMl) as 

= Jdx 1 f dx 2 J dp^ip^if^f^^)] 

x^k? 2 -(P«+P*) 2 ) (66) 
q y 
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Integrating over q on both sides of the above equation, we obtain 

d pd {p,p d ) r. q/ , r q ( .Ana 2 2 

= ^) r ^)]^rt, (67) 

where the q contribution is isolated in 

rLiv(* 2 ) = / dp NPpd {p N )[^fl N {^)\ . (68) 
J Pn Pn 

The derivation of Eq.( l68l can be extended to have q in deuteron and q in the projectile 
proton. We finally obtain 

= Y.if q M)\ E n d A^)\ + f«(xi)[ E Wl^-j • (69) 

12 (jr jV=p,Tl jV=p,7l H 

We use the charge symmetry to calculate the parton distribution functions for the neutron 
from that of proton: f£ = fp, f% = fp, ft = fp , fn = fd- Furthermore P PcL (pn) is the same 
for neutron and proton. Including the charges for u and d quarks appropriately, Eq. flB^]) can 
be written as 



D. Calculation of p Pd {pN) 

To calculate the p-d DY cross section Eq.( )42l) . we need to generate the nucleon momentum 
distribution pp d (pW) defined by Eq.(l39l). Since the available realistic nucleon-nucleon poten- 
tials only provide the wavefunctions in the rest frame, we boost these wavefunctions within 
the relativistic quantum mechanics proposed by Dirac, as reviewed in Ref.(32[. We note here 
that the parton distributions defined by Eqs. (122]) - (f23|) or Eqs. fl56|) -f l57|) are also within the 
same framework. We thus have a consistent relativistic description of both the nucleonic 
structure of the deuteron and the partonic structure of the nucleon. Such a consistency is 
essential in investigating the Fermi motion effect on DIS and DY processes. 

Starting with a mass operator (Hamiltonian in the rest frame) which defines the mass and 
the wavefunction Xd{k) of the deuteron in its rest frame, one can generate the wavefunction 



$ Pd of a moving deuteron with momentum p^ using one of the three possible forms [3 ll. |32| of 
the relativistic quantum mechanics. Furthermore, one can identify Xd{k) with the deuteron 
wavefunctions generated from the available realistic nucleon-nucleon potentials, as discussed 



in Refs.|33l. l34j. Here we consider the Instant-Form (IF) and Light-Form (LF). The formula 
for calculating p Pc1 (pn) can be derived by using the procedures presented in Ref.{32[ and will 
not be detailed here. We will simply write down the formula used in our calculations. 



1. Instant Form 



The deuteron wavefunction Xd(k) in the pd = rest frame is normalized as 

dk\ Xd (k)\ 2 = l. 
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If we neglect the interaction in the Lorentz Boost and follow the procedures detailed in 
Ref. 32], it is straightforward to show that Eqs. flM]) and (I39I) lead to 



PpM 



< %d 








Pd 


^Pd ^ 



d(p,p 2 ) 



X d (k)\\ (71) 



with 

I d(p d , k) _ uj m (k)uj m (k)u Mo (p d ) 
d(p,p 2 ) u m (p)u m (p 2 )M 



(72) 
(73) 

where m is the nucleon mass, u; m (k) = \J m 2 + k 2 , and 

P2 = Pd-P, (74) 

% = p+ w [ w& - Umm ' (75) 

H = u m {p) + u m (p 2 ) , (76) 
Mo = [H 2 - P 2 d ] 1/2 . (77) 

2. Light Form 

The Light-Form (LF) momentum is defined by 

P = (P + ,P±) , 

where 

p + = u m (p) +p 3 , 
P± = (p\p 2 )- 

By using the above definitions for all momenta in Eqs. fl34l) and (1391) and following the 
procedures detailed in Ref. 3^], one can show that the nucleon momentum distribution in 
LF is 



< ®pM b p 


^Pd 


- ^p d 


^p d ^ 



p L p?(p) 

,d(p d ,k) T 2 



d(p,p 2 ) 



Xd(k)\\ (78) 



where 

d(p d ,k) _ u m (k)u m (k)p d ~ 



(79) 

d(p,p 2 y p + ptMo 
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FIG. 3: Nucleon momentum distribution p Pd (p) of a deuteron moving with a momentum = |p^| 
in p-cZ collisions at the proton laboratory energies (solid ), 120 (dashed), 800 (dash-dotted) 
GeV 



with 



P2 


= Pd-P, 






k± 


pi 

= —p± - 

Pd 


P + -, 

- —P2± , 
Pd 


k 3 


= -(p + - 


pi)\ 


m 2 + kj_ 




M 


= 2u m (k) 







To use LF results in the calculation of DY cross section Eq. 
the LF variables into the usual variables. We thus have 



PpM 



P 



OJr. 



(80) 
(81) 

(82) 
(83) 

we need to transform 
(84) 



VI. NUMERICAL RESULTS 

To calculate p+d DY cross sections, we need to first determine the average nucleon 
momentum p z ave in Eq.f l68|) . This is done by examining the dependence of the nucleon 
momentum distribution p Pd (p) on the deuteron pa using the formula given in subsection 
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V.D. In Fig|3j we show p Pd (p z ) = I Pp d {p)dp x dp y of the Instant Form (IF) which are used 
in our calculations for the p-d DY cross sections at E = 0, 120, 800 GeV. The Argonne-V18 
potential |42j is used here to generate the deuteron wavefunction Xd(k)- We see that as 
Pd increases, the spreading of the momentum distribution is wider. This is of course the 
consequence of Lorentz contraction of the density distribution in coordinate space. In the 
same figure, we also indicate the deuteron momentum pd for each case. We observed that 
their peaks are near pa/2 at each energy, and we hence set p z ave = p d /2 in the calculation of 
DY cross sections using Eqs. fl68|) and flTOj) . 

We have also carried out calculations for the Light-Form (LF) and find that the results 
are indistinguishable from that shown in FigJSJ This turns out to be due to the fact that in 

the region where |x(^)| 2 dominant, the Jacobin of IF and if^^l of LF are almost 

the same numerically. 

We use the CETEQ5m[4l[ parton distributions to perform calculations using Eq.( l60l for 
p-p and Eq.f l6"£j) and Eq. fjTUj) for p-d DY cross sections. If we neglect the p^-dependence 

and set p z N — >• p ave such that 



-> ftrfa) since J dp N p Pd (p N ) = 1 as 



in the parton distribution function f q (x2P ave /p z N ) in Eq. 
f 9 ( x 2Pave/PN) ~> Pfa), we then have T q pdN (x 2 
defined by the normalization of deuteron wavefunction Xd- The calculations based on this 
simplification are then identical to those from using Eqs. of Ref. Q . Clearly the results 
from using this frozen-nucleon approximation do not include the Fermi Motion (FM) effects. 
Our interest is to see how these "No FM" results are compared with the "With FM" results 
calculated from using the convolution formula Eqs. (1681) and (ITOjl . 
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FIG. 4: Ratio a pd /a(pp) at E = 800GeV. The data are from Ref.|j|. 



We first investigate the data at 800 GeV of Ref. (B] by calculating the ratio 



R 



pd/pp 



dxidx 



p-) /(2x *np 



dx\d%2 



(85) 
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FIG. 5: Ratio a pd /a(pp) at E = 120GeV and x x = 0.6. 



Our results are shown in FigJH We see that the Fermi motion has less than 1% effect on 
the ratio a pd /(2a pp ) in the low x 2 < 0.25 region covered by this experiment. Our "No FM" 
results are similar to that presented in Ref.[5j and agree with the data. 

In FigfSl we show the Fermi motion effect at 120 GeV which will be considered in the 
upcoming experiment at Fermi Laboratory. The results are for x\ = 0.6. We see that at 
small x 2 , the FM effect is also very small as in FigJUfor 800 GeV. However the FM effects 
become increasingly larger with x 2 \ about 25% as x 2 — > 0.6. This is similar to what is known 
in the DIS studies. 

To see the Fermi motion effects more clearly, we present results for x\ = 0.9 in Figj6j 
We see that the ratio a pd /(2a pp ) — > about 2.5 as x 2 —> 0.9. This can be understood from 
FigJT] where we show the differences between the convoluted parton distribution rjj ^(^2) , 
defined by Eqs. (l6"8~|) . and /^(a^) of CETEQ5m. We see that their differences can be as large 
as a factor of about 2 as x 2 — > 0.9. Similar differences are also for other parton distributions. 

The results presented in FigsfSE are from using the deuteron wavefuncition x(k) gener- 
ated from Argonne-V18 NN potential. Similar results have also been obtained by using the 
wavefunctions from Nijmeganj43j, CD-Bonn 44], and Paris 45] NN potentials. 



VII. SUMMARY AND DISCUSSIONS 



Within a formulation based on relativistic quantum mechanics, we have developed con- 
volution formula for calculating the p-d Drell-Yan cross sections using the available parton 
distribution functions |4l| and the nucleon momentum distributions generated from realistic 

When the Fermi motion in the deuteron is neglected, our 



nucleon-nucleon potentials[42-45 



formula are reduced to the factorized form Eq.(2) which was used in the previous analysis. 
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FIG. 6: Ratio a pd /a(pp) at E = 120GeFand x x = 0.9. 



By comparing the results from the convoluted and factorized calculations, we have shown 
that in the small x 2 < 0.3 region, the Fermi motion effect is very small and our results for 
the ratios R pd/pp = a pd /(2a pp ) agree well with the data at 800 GeV of Ref . (BJ . On the other 
hand, the Fermi motion effect can influence significantly the p-d Drell-Yan cross sections in 
the large Bjorken x 2 > about 0.4 region. At 120 GeV the predicted Fermi Motion effect can 
enhance the ratios R p d/ PP by about 20 % at x 2 ~ 0.6 and about a factor of 2.5 at x 2 — >■ 1.0. 
Our results suggest that the Fermi motion effect must be included in extracting the d/u 
ratio in the nucleon from the upcoming experiment at Fermi Laboratory. 

In this work, we have only investigated the most trivial nuclear effects on DY processes. 
To extract the d/u ratios precisely, we need to extend our formulation to account for other 
nuclear effects. By using the approximations Eqs. ([2"Tj) and ( 12"5"|) for the time- components 
of 5-functions in our derivations, the binding effects on the partons in the nucleon and the 
undetected fragments X p and X N are not treated rigorously. To improve this incomplete 
treatment of binding effects, one must depart from the simple convolution formula by per- 
forming calculations using the spectral functions of the nucleon in terms of parton degrees 
of freedom. Such information is not available at the present time. 

The use of impulse approximation, as illustrated in FigfS] , only accounts for the leading- 
order term of the full DY amplitude. It is also necessary to investigate the accuracy of 
this approximation. For example, we perhaps need to investigate the double scattering 
mechanisms which are considered in Ref. |48| to describe the shadowing effects on DIS. We 



also need to investigate the effects due to pion excess, as discussed in Refs. |25|, [50 



Finally, we mention that there exists other relativistic formulations of DIS and DY pro- 
cesses. The differences between different approaches are mainly in describing the structure 
of the considered composite systems, nucleons or nuclei, in terms of interactions between 
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FIG. 7: Convoluted and factorized u distributions at E = 120GeV and x\ = 0.9. The results for 
d are very similar and are therefore not shown. 



their constituents. In the formulation considered in this work, all particles are on their 
mass shells, but the systems can be off the energy shell during interactions. Thus it is also 
necessary to examine how our formulation differs from other formulations within which the 
nucleons in the deuteron can be off their mass-shell p 2 N ^ m%, such as the covariant formu- 
lation given in Ref . |49| . Similar theoretical questions also exist in many-years investigations 
of nuclear dynamics at intermediate and high energies; such as the differences 5l(] between 
different nucleon-nucleon potentials constructed using different three-dimensional reductions 
of Bethe-Salpeter equations. 
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